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We evaluate the (9(^ HL ) correction to the rate of B — > X s 7 decay, i.e. we resum all the 
^ ! (D[(a em \n Myy / mf) x (a s lnM^/mj)"] corrections for n = 0, 1,2, ... . Our calculation differs 
^ | from the previously available one by that it takes into account the complete relevant set of 
operators. The correction is found to be negligible, i.e. it is below 1%, in accordance with 
the former results. 



Abstract 



* This paper is based on the M.Sc. thesis of the first author. 



At present, the next-to- leading logarithmic QCD analysis of B — > X s "y decay allows 
predicting its branching ratio with around 10% accuracy 0. The current accuracy on the 
experimental side is around 15% [[|, which is expected to be improved soon. Therefore, it is 
important to examine the size of the dominant electroweak corrections to this decay mode. 

Prior to an explicit calculation, two sets of electroweak corrections are expected to be 
dominant. The first set consists of corrections that are enhanced by the ratio rn^/M^. Such 
contributions have been evaluated in ref. || and found to have less than 1% effect on the 
branching ratio, owing to an accidental cancellation. The second set consists of corrections 
that are enhanced by the large logarithm In M^/rn^. The authors of ref. [Q have calculated 
them and found a contribution of only around 1% to the branching ratio. 

However, when the O (a em In M^/mf) correction is small, it might happen that some of 
the 0[(a em In j mf) x (a s In / 7n|) n ] corrections are much bigger, while only the n = 
case was included in ref. f|] . Since a s In / m\ is close to unity, a naive estimate for the 
sum of all the 0[(a em \ia M^r / mf) x ( a s h^M^/m^)™] corrections is - x (number of 
order unity) = 6.5% x (number of order unity). The above estimate is close to the overall 
0(10%) uncertainty of the prediction for the branching ratio. Thus, there is a good reason 
for performing a complete 0( 9 ^ B -) calculation. 

Resummation of the large logarithms In / m\ from all orders of the perturbation series 
is most conveniently performed in the framework of an effective theory that is obtained from 
the Standard Model by decoupling the top quark and the heavy electroweak bosons. 

The part of the effective theory lagrangian that is relevant to b — > 57 reads 



icy 

d,s,c,b) + -^-V t * s V tb 



A=l i=3 

where are elements of the CKM matrix, while Cj(/i) and Cf([i) are the Wilson coefficients 
at the following operators: 



Pi = 


(s L ^T a c L )(c L rT a b L ), 


P2 = 


(sLlt,c L )(c L Yb L ), 


Ps = 


(sLi^b^Eqiqi^q), 


Pa = 


(s LllM T a b L )J2 q (qYT a q), 


P 5 = 


{sLl^lM £ ? (?7 M W 3 <z), 


Pe = 


(si7 w 7M 2 7M3^ a ^)E g (97 m T M2 T M3 ^°g), 


pQ _ 


(sz7A)£ 3 Qq(q^q)> 


pQ _ 


(s Ll ,T a b L )E q Q q (qrT a q), 
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if = (WnWM E 9 Q q (q7^ 2 r 3 q), 

P? = (s L ^^ :i T a b L ) E q Q q (qY^r s T a q), 

Pr= ^m b {s L <7^b R )F^, 

Ps= j^m b (s L a^T% R )G%. (2) 

The small CKM matrix element V ub as well as the s-quark mass are neglected here. 

The above set of operators closes under the QCD and QED renormalizations.0 It is the 
QED renormalization that forces us to introduce the operators P®, in which sums over 
flavours are weighted by their electric charges. These operators were absent in the QCD 
analysis of ref. JL], because O^ 9 ^) corrections were neglected there. 

In order to evaluate the Oi^^) terms, one has to calculate the anomalous dimension 
matrix for all the above operators up to order 0(a em ), and then solve the Renormaliza- 
tion Group Equations (RGEs). Such a calculation has already been performed in ref. |5|]. 
However, the operator basis was truncated there to {Pi, P2, P7, Ps} only. 

The main purpose of the present paper is a complete evaluation of the 0{ g ^ Sk ) corrections, 
i.e. including all the relevant operators from eq. (0). At the same time, we shall check the 
results of ref. || and verify whether truncating the operator basis was a good approximation 
there. 

Evaluating the anomalous dimension matrix at 0(a s ) and 0(a em ) proceeds in full analogy 
to the well-known calculations of the leading-logarithmic QCD effects in B — > X s 7 (see 
[0 and references therein). One needs to find the one- loop mixing among the four-quark 
operators, the one-loop mixing in the {P7, Pg} sector, as well as the two-loop mixing of the 
four-quark operators into P7 and P§. Diagrams with virtual photons need to be included as 
well. 

As in the case of the former QCD analyses || [I], it is convenient to introduce the so-called 
"effective coefficients" before the RGEs are solved. They are given by the following linear 
combinations of the original Wilson coefficients: 



A=3 *i 



a(/i) - |c?(a*) 

C t (p) - §Cf (//) 



for i = 7, 

for i = 8, (3) 



Ci(ji), otherwise. 
The numbers and zi are defined so that the leading and the Oi^^) contributions to the 

b — > 57 and b —> s gluon matrix elements of the effective hamiltonian are proportional to 

2 It closes off-shell, up to non-physical operators that either vanish in four dimensions or vanish by the 
QCD x QED equations of motion. The existence of leptons is ignored here, because their effect b — ► sj is 
of higher order in QED. 
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the corresponding terms in Cj^ and Cf^ , respectively. In dimensional regularization with 
fully anticommuting 75, we have y = (0, 0, — |, — |, — y , — tt) and 2 = (0, 0, 1, 



20, -m. 



The effective coefficients evolve according to their RGEfJ 
driven by the anomalous dimension matrix 7 e -^(/i) 



a—C e - 
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4vr rs 4vr rem 



(4) 



(5) 



The matrices ^^ff and 7^ e ^ are regularization- and renormalization-scheme independent, 
contrary to the matrices governing the evolution of the original coefficients Ci{fj). 
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(O)eff 



*(0)eff 
1 em 



s for 


7 (0)e// 


and 


A(0)e// 
/ em 


are 


the following 












Pi 


P'2 


Pa 


Pi 


Ps 


P 6 


tdQ 

r % 


pQ 


tdQ 
^5 


pQ 
^6 


Pi 


Ps 


-4 


3 





2 
9 




















208 
243 


173 

162 


12 








4 
3 




















416 
81 


70 
27 











52 
3 





2 














176 
81 


14 

27 





- 


40 

9 


100 
9 


i 

9 


5 

6 














152 
243 


587 
162 











256 
3 





20 














6272 
81 


6596 
27 





- 


256 
9 


56 
9 


10 
9 


2 
3 














4624 
243 


4772 
81 











8 
9 











-20 





2 


176 
243 


14 
81 











16 








_40 


_52 


4 


5 


136 


295 


27 


9 


3 


9 


(5 


729 


486 











128 
9 











-128 





20 


6272 
243 


764 
81 











184 








256 


160 


40 


2 


39152 


1892 


27 


9 


3 


9 


3 


729 


243 
































32 
3 



































32 
9 


28 
3 


8 
3 

















32 
27 











832 
729 


22 
243 





8 
3 














8 
9 











208 
243 


116 
81 




















76 
9 





2 
3 





20 

243 


20 
81 




















32 
27 


20 
3 





2 
3 


176 
729 


14 
243 




















496 
9 





20 
3 





22712 
243 


1328 
81 




















512 
27 


128 
3 





20 
3 


6272 
729 


1180 
243 
































16 
9 


8 
3 



































8 
9 



Pi 

P-2 
P 3 
Pi 
P, 

Pe, 
P Q 

-<3 

pQ 
pQ 
P Q 

Pi 
Ps 

Pi 

Pi 
P 3 
Pi 
P, 
Pe 
Pi 
Ps 



3 In eqs. (Q)— (^), we do not distinguish between Ci and Cf, i.e. we write these equations as if the 
operators were numbered from 1 to 12. 
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Details of their evaluation can be found in ref. |7|]. 

The rows corresponding to P®,...,P® in the matrix 7^ e ^ have not been given explicitly 
above. They would be relevant only at higher orders in a em , because the coefficients of these 
operators start at O(^). 

In the above matrices, the QED mixing of P 3 ,...,P 6 into P7 and Pg as well as the QCD 
mixing of P3 ,...,Pg into P7 and Pg are given for the first time. As far as the remaining 
entries are concerned, our results agree with the old ones of refs. [§, [| and ||. However, in 
order to perform a comparison, one needs to make a linear transformation of our matrices 
to the "old" basis of the four-quark operators used in those articles. 

The solution of the RGE (£|) with initial conditions at /i = M\y has the following form: 

cf f {n) = c^ ff M + ^cr (0)e// (/i) + ... , (6) 

a a {p) 

where 

Cf )eff {n) = Y,Vikri ak V k fCp eff (M w ), (7) 
Cr i0)eff M = J_ E fjf k ~\ k (V-\lSl eff ) T v) kl (8) 

^0 k,l,j ™ ' 



(3q = y and 77 = . The matrix V" and the numbers a k are obtained via diagonalization 



f and - - 
of (7i 0)e// ) T 

(v-\t )eff ) T v) kl = 2(3 a k 5 kl . (9) 

The above solution to the RGE is identical to the one found in refs. || and 0.0 

Note that the ) terms in the Wilson coefficients vanish at the initial scale \x = M^y. 

It must be so, because no such terms can arise from perturbative matching of the SM and 
the effective theory amplitudes. The relevant matching conditions are thus the same as the 
leading-order ones in ref. JO: 



cP eff (M w ) = C ( ?\M W ) = < 

C? (0)e// (M W ) = 0, 

where x = mj/M^. 

4 Except for the misprint r]i «-► rjj in eq. (A. 10) of ref. || 



0, for i = 1,3,4,5,6, 

1, for i = 2, 

3x 3 -2x* i 1 -8x 3 -5x 2 +7x for i = 7 (10) 



-3x- 2 !„ „ i -x 3 +5a: 2 +2x 



4(a _^lnx+ ^ 8 y_y , for 2 = 8 
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After substituting the explicit anomalous dimension matrices and the initial conditions 
to the solution of the RGE, we find 



Cf e// (A*) 

C (0)e// (^) 



i=i 



( 88 i6 40 7 32 9 \ -fou,, x 



+ 



32 



704 16 640 14 32 _ x 

77 23 -| 77 23 H 77 23 

1725 ' 1449 ' 1449 ' 575 



526074716 _47 65590 _m 

77 23 -| T] 23 

4417066408125 ' . ' 



V -^ci 0) (M w ) 
1686113 ' -t (l>'>»" ■>>>»■' ' 



i=i 



'12) 



with the values of a;, hi, h\ and ^ given in table 1. 



i 


Oj 


hi 


K 


K 


1 


14 

23 


626126 
272277 


50090080 
131509791 


10974039505456 
21104973066375 


2 


16 

23 


56281 
51730 


107668 
646625 


13056852574 
29922509799 


3 


6 

23 


3 
7 


3254504085930274 
23167509579260865 


718812 
6954395 


4 


12 

23 


1 

14 


34705151 
143124975 


154428730 
12196819523 


5 


0.4086 


-0.6494 


-0.2502 


-0.1374 


6 


-0.4230 


-0.0380 


0.1063 


-0.0078 


7 


-0.8994 


-0.0186 


-0.0525 


-0.0023 


8 


0.1456 


-0.0057 


0.0213 


-0.0001 



Table I. The numbers dj, hi, h\ and h- entering eqs. ( pTTj) and ([12]). 
Setting fib to 5 GeV and taking the remaining parameters from ref. 0, one finds 77 ~ 0.56, 



mj/M^r ~ 4.7 and a em / a s (fib) — 0.036, which implies 



c e 7 



(0)eff 
7 

:m (0)e//^ 



and, in consequence, 



ABR[B -> X sl ] _ 2 _ a em 



= -0.313, 
= 0.033, 

cv em(0)e// (/i b ) 



(13) 
(14) 



-0.8%. 



(15) 



BR[B^X sl ] ~aM cf e// (/i b ) 
The above 0(^ m -) correction to BR[B — > Xgj] appears to be much smaller than the simple 
estimate presented at the beginning of this article. The structure of eqs. ([Ll]) and (li~2f) leads 
to a naive expectation that cf^ e ^ (fib) and C^ m ^ e ^ (fib) are similar in magnitude. However, 
after the numerical evaluation, one of them turns out to be almost ten times smaller. 
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Although our analytical result for (jvm,(o)ef} (fib) is different from eq. (11) in ref. 0, the two 
formulae are numerically very close, for realistic values of rj. Thus, truncating the operator 
basis is a correct approximation in the present case. 

To conclude: We have performed a complete calculation of the 0( 9 ^ !L ) correction to the 
branching ratio of B — > X s 7. The correction is found to be approximately equal to —0.8%, 
i.e. negligibly small. However, without an explicit calculation, a correction several times 
larger could not have been excluded. 
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